Transport properties of strongly correlated electrons in quantum dots using a simple 

circuit model 
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Numerical calculations are shown to reproduce the main results of recent experiments involving 
nonlocal spin control in nanostructures (N. J. Craig et al., Science 304, 565 (2004)). In particular, 
the splitting of the zero-bias-peak discovered experimentally is clearly observed in our studies. To 
understand these results, a simple "circuit model" is introduced and shown to provide a good 
qualitative description of the experiments. The main idea is that the splitting originates in a Fano 
anti-resonance, which is caused by having one quantum dot side-connected in relation to the current's 
path. This scenario provides an explanation of Craig et al.'s results that is alternative to the RKKY 
proposal, which is here also addressed. 

PACS numbers: 71.27.+a,73.23.Hk,73.63.Kv 



The observation of the Kondo effect in a single quan- 
tum dot (QD)i and the subsequent theoretical and ex- 
perimental studies of more complex structures, such as 
two QDs directly coupled through a tunable potential 
barrier—, has provided impetus for the analysis of more 
elaborate systems. In a recent seminal work, Craig et 
al£ report on the possible laboratory realization of the 
two-impurity Kondo system. Two similar QDs are cou- 
pled through an open conducting central region (CR). 
A finite bias is applied to one of the QDs (QD1 from 
now on) as well as to the CR, while the other QD (QD2) 
is kept at constant gate potential. The differential con- 
ductance of QD1 is then measured for different charge 
states of QD2 and different values of its coupling to the 
CR. The main result was the suppression and splitting 
of the zero-bias-anomaly (ZBA) in QD1 by changing the 
occupancy of QD2 from even to odd number of electrons 
and by increasing its coupling to the CR. A Ruderman- 
Kittel-Kasuya-Yosida (RKKY) interaction between the 
QDs was suggested as an explanation for the observed 
effects^. The importance of Craig et al.'s experiments 
cannot be overstated: the possibility of performing nonlo- 
cal spin control in a system with two lateral QDs has po- 
tential applications in QD-based quantum computing 5 . 

In this Letter, numerical simulations in good agree- 
ment with the experiments are presented. The central 
conclusion of this work is that our computational data, 
and as a consequence the experimental results, can be 
explained using a very simple "circuit model" , where one 
of the elements is a T-connected QD that has an intrin- 
sic reduction of conductance with varying biases. This 
proposal is an alternative to the more standard RKKY 
ideas. Fig. la depicts the experimental set up used in 
the measurements of Craig et al? with the labeling used 
in this Letter. Figure lb is a schematic representation 
of the system, introducing two different tunneling pa- 
rameters (hopping matrix elements t' and t") and the 
Coulomb repulsion U in each QD (assumed the same for 
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FIG. 1: (a) Experimental setup used in Ref. |3j (b) Illustra- 
tion of the model studied in this Letter (see text for details). 



simplicity). To model this system, the Anderson impu- 
rity Hamiltonian is used for both QDs: 



H d = 
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where the first term represents the usual Coulomb repul- 
sion between two electrons in the same QD, and the sec- 
ond term is the effect of the gate potential V g i over each 
QD. QD1 is directly connected to the left lead and to the 
CR with hopping amplitude t', while QD2 is connected 
only to the CR (with hopping amplitude t"), which itself 
is connected to the right lead with hopping amplitude t 
(which is also the hopping amplitude in both leads, and 
our energy scale). In summary, 
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where c\ ia (cl ia ) creates an electron at site i with spin a 
in the left (right) lead. The CR is composed of one tight- 
binding sitefi, unless otherwise stated. Site '0' is the first 
site at the left (right) of QD1 (CR) in the left (right) lead. 
The total Hamiltonian is Ht = + -ffi ca ds + H12 ■ Note 
that for V g i = V g 2 = —U/2, the Hamiltonian is particle- 
hole symmetric. To calculate the conductance G, using 
the Keldysh formalism 7 , a cluster containing the inter- 
acting dots and a few sites of the leads is solved exactly, 
the Green functions are calculated, and the leads are in- 
corporated through a Dyson Equation embedding pro- 
cedure. Details of the embedding have been extensively 
discussed before^. All the results shown were obtained 
for U — 0.5, t' = 0.2, zero-bias, and zero temperature. 




FIG. 2: (a) Variation of G with V s i in QD1 (solid curves) 
and variation of (712) (occupancy of QD2 per spin orientation 
- dashed curves) for t" = 0.2 and three different values of 
V g 2. For I42 = —2.0 (red), QD2 is occupied by 2 electrons 
(("■2) = 1) for any Vgi and the conductance through QD1 is 
essentially the same as if QD2 was not present. For higher 
values of V g 2, the average value of (712) decreases and becomes 
dependent on V g i (decreasing for lower values of Vgi). This is 
accompanied by a suppression of the ZBA (for —0.35 (green)) 
and also by a splitting of the ZBA (for V S 2 = —0.3 (blue)), 
(b) Variation of G and (n 2 ) with t" (0.2, 0.4, 0.5) at a fixed 
value of Vg2 = —0.5. As the value of t" increases, the average 
value of (712) decreases and this is again accompanied by a 
suppression of the ZBA. (c) Same as in (b), but now for V g 2 = 
—0.25 (particle-hole symmetric point) and t" = 0.0, 0.1, and 
0.2. Note that G vanishes at V g i = —0.25, where (712) = 0.5, 
for all finite values of t". 

In Fig. 2, results for the conductance across QD1 (solid 
curves) and for the occupancy per spin orientation (712) of 
QD2 (dashed curves) are presented. In Fig. 2a, t" = 0.2 
and V g 2 varies from —2.0 to —0.3. For V g 2 = —2.0 (red), 
QD2 is occupied by two electrons ((712) = 1) and the 
conductance of QD1 displays the characteristic Kondo 
behavior reported before^. For V g 2 = —0.35 (green) the 
average value of (712) decreases to « 0.7 (« 1.4 electrons 
in QD2) and (712) now depends on V g \. In addition, G 



decreases in comparison to the result obtained for V g 2 = 
—2.0. Then, these numerical results are qualitatively in 
agreement with the experimental results shown in Fig. 2 
of Craig et al£, namely, by decreasing the occupancy of 
QD2, from even to odd number of electrons, the ZBA in 
QD1 is suppressed. As V g 2 is further increased (—0.3 
(blue)) a qualitative change occurs: For values of V g \ 
where (712) ~ 0.5 (QD2 singly occupied), the conductance 
of QD1 vanishes and therefore there is a narrow dip in 
G. This splitting of the ZBA is remarkably similar to 
that observed in Fig. 3A of the experimental results^. 
For finite-temperature calculations, the dip in G will not 
reach zero, resembling even better the experiments^. 

To further test the similarities between simulations and 
experiments, in Figs. 2b and 2c results for G and (712) 
are shown for fixed V g 2 and different t" values. In Fig. 
2b, where V g 2 — —0.5, as t" increases from 0.2 to 0.5 
there is only a slight decrease of G. This is accompanied 
by a slight decrease in the average value of (712), from 
ss 0.9 to « 0.7. A more dramatic change is obtained 
in Fig. 2c, where V g 2 = —0.25, and t" varies from 0.0 
to 0.2. By increasing t" from 0.0 (red curves) to 0.1 
(green), the ZBA is now split in two and (712) acquires a 
dependence on V g \. As t" further increases (0.2 (blue)), 
the dip becomes wider, the two side-peaks decrease and 
G still vanishes for (77,2} = 0.5 (one electron in QD2). 
Our calculations show that, if (77,2) varies around 0.5, 
the dip in G is present for all finite values of t", with 
a width proportional to t". Comparing the results in 
Figs. 3A and 3B of Craig et al£ with Figs. 2c and 2b in 
this Letter, respectively, one notices a striking similarity: 
The splitting of the ZBA observed in the experimental 
results (their Fig. 3A) , when the number of electrons in 
the control QD is odd and the coupling to the central 
region is increased, is very similar to the dip in G for 
all finite-i" curves in Fig. 2c (as mentioned above, at 
finite temperatures, one expects that the dip in G will 
not reach zero). When the occupancy of QD2 is even 
(Fig. 3B in the experimental results^ and Fig. 2b in this 
Letter), the G dependence on t" is much less significant 
and the splitting of the ZBA does not occur. 

What is the origin of these results? Below, it will be 
argued that a qualitative description of the results can be 
achieved by analyzing the two quantum dots through a 
so-called 'circuit model'. This model starts with the con- 
ductance of each QD calculated separately, as indepen- 
dent elements of a circuit, and then the conductance of 
the 'complete circuit' is obtained by combining the con- 
ductances of the two elements connected in series. Fig. 
3 describes schematically the steps involved in this ap- 
proach. In Fig. 3a, the complete system formed by QD1 
and QD2 (shown in Fig. lb) is divided into two compo- 
nents. QD1 is modeled as a QD connected directly to left 
(L) and right (R) leads, while QD2 is modeled as a side- 
connected QD 12 . Fig. 3b shows the respective conduc- 
tances and occupancies for each independent element vs. 
gate voltage, and Fig. 3c represents the scattering pro- 
cesses (represented by transmission and reflection ampli- 
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FIG. 3: Schematic representation of the main ideas behind 
the "circuit model". In (a), the system represented in Fig. 
lb is divided into its constituent elements: QD1 is modeled 
as a QD connected in series with the leads and QD2 is mod- 
eled as a side-connected QD. The curves in (b) represent the 
conductance and occupancy of each separate circuit element 
vs. the applied gate potential, (c) Schematic representation 
of how the two individual elements are connected to form 
the final 'circuit': Incident and reflected wave amplitudes are 
represented in the right side of QD2 by black arrows. A trans- 
mitted wave through QD2 undergoes multiple reflections be- 
tween the two quantum dots until it is finally transmitted 
past QD1. The superposition of all these processes results in 
the final conductance for the 'circuit'. 



tudes) that an electron undergoes while moving through 
the complete 'circuit'. The superposition of all these pro- 
cesses leads to the total transmittance (proportional to 
the conductance) for the circuit model. This can be cal- 
culated in two ways: coherently or incoherently^. Since 
there is no qualitative difference between them, and in 
order to keep the simplicity of the model, we present the 
incoherent results. The equation which provides the final 
transmittance for the processes depicted in Fig. 3c is 



T = 



TiT 2 



(4) 



where the transmittances T\ and T 2 are proportional to 
the conductances for QD1 and QD2, as depicted in Fig. 
3b, andi?i(2) = 1— 2\(2) are the reflectances. To calculate 
T, one needs to establish how T 2 depends on Vgi. The 
natural way to do that is to use the dependence of (712) 
on V g i, as depicted in Fig. 2, and then use the relation 
between conductance and occupancy, as shown in the red 
curves in Fig. 3b. In other words, the functional relation 
can be expressed as T 2 = T 2 ((n 2 )(V g i)). It is not sur- 
prising that in a strongly correlated system like the one 
being analyzed here, the variation of the gate potential 
of QD1 will influence the charge occupancy of QD2, and 
in turn this will influence the conductance through QD1. 

In Fig. 4, conductance results using Eq.(4) are shown 
for the same parameters as in Fig. 2. Although the quan- 
titative agreement varies, there is good overall qualitative 
agreement. All the trends are correctly reproduced and 
some of the details are quite similar, such as for exam- 
ple the asymmetric shape of the curves at higher values 
of V g2 (—0.35 and —0.3) in Fig. 4a. It is important to 
notice that there are no adjustable parameters in the cir- 
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FIG. 4: Same as in Fig. 2, but now using the 'circuit model' 
for the calculations. 



cuit model here presented. The only input necessary is 
(n 2 ) vs. Vgi, which is obtained through a calculation for 
the complete system (values displayed for (n 2 ) in Fig. 
2). The success of the circuit model implies that the dip 
in G arises from the Fano anti-resonance which cancels 
the conductance of QD2 (red solid curve in Fig. 3b). 
The Fano anti-resonance can be seen as a destructive in- 
terference process between two different trajectories an 
electron can take on its way to QD1: it can cross the 
CR without passing through QD2; or it can visit QD2, 
return to the CR and then proceed to QD1— . 

The similarities between the experimental results and 
our simulations suggest that our model and numerical 
technique have captured the essential physics of the ex- 
periments. However, these same experiments have also 
been explained using RKKY ideasi Can our numerical 
results be also understood in this alternative context? 
To try to answer this question, several calculations were 
performed with different parameter values and number 
of sites in the CR 14 . In Fig. 5a, results for spin cor- 
relations between QD1 and QD2 (denoted Si • S2) are 
presented for the same parameters used in Fig. 2c. At 
t" = 0.0 (red curve) QD1 and QD2 are uncorrelated as 
expected. As t" increases to 0.1 (green), and then 0.2 
(blue) , it is observed that in the region where G reaches 
its maximum value (see Fig. 2c), Si • S2 also assumes a 
maximum value and it is positive (ferromagnetic (FM)). 
For t" > 0.2 (not shown), Si • S2 saturates and starts 
decreasing. The maximum of Si • S2, for all values of t", 
decreases even further as the size of the central region in- 
creases (the results in Fig. 5a are for a CR with just one 
site). In addition, the sign of Si ■ S 2 alternates as the 
size of the CR increases and the QDs move farther apart 
from each other. In Fig. 5b, results for the spin corre- 
lation between QD1 and its neighboring site in the CR 
(denoted Si • S c ) is shown for the same parameters as in 
Fig. 5a. Si • S c is a rough measure of the Kondo corre- 
lation in QD1, having a direct connection with the ZBA 
in Fig. 2c. Indeed, for t" = 0.0 (red) when G reaches 
the unitary limit, a robust antiferromagnetic (AF) cor- 
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relation develops between QD1 and its neighboring site 
in the CR. For t" = 0.1 (green), despite the narrow dip 
in G, the side-peaks are still close to the unitary limit 
(see Fig. 2c) and Si • S c is still strongly AF. However, 
for t" = 0.2 (blue), both G and Si • S c are strongly sup- 
pressed, in qualitative agreement with a suppressed ZBA 
due to a weakened Kondo resonance. 
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FIG. 5: (a) Spin correlation Si ■ S2 between QD1 and QD2 
for the same parameters as in Fig. 2c. For t" = 0.0 (red 
curve), the two QDs are uncorrelated (Si ■ S2 =0). For finite 
t" (0.1 (green) and 0.2 (blue)), Si ■ S 2 is FM and reaches 
its maximum value in the region where G is maximum, (b) 
Kondo correlations Si • S c between QD1 and the central site 
for the same parameters as in (a). All values are AF and 
they decrease in amplitude as t increases, underscoring the 
decrease of the Kondo effect as the FM correlation between 
QD1 and QD2 increases (compare with (a)), (c) Variation of 
G as Ui (Hubbard interaction in QD1) assumes the values 0.4, 
0.2 and 0.0. Note that the dip in G becomes slightly narrower 
as Ui decreases, however it does not disappear. 

The results thus far seem to indicate that the CR could 
be mediating a long range coupling between QD1 and 
QD2, with the characteristics of an RKKY interaction. 
However, the magnitude of the maximum value of Si • S2 
(see scale in Figs. 5a-b) is too small to account for all the 
effects observed in the conductance in Fig. 2c. One pos- 
sible way of increasing Si ■ S2 is by coupling QD1 more 
strongly to the CR than to the left lead. This was ex- 
actly the setup chosen in Ref. 0, where those authors 
performed the measurements with asymmetric couplings 



to the left (Fl) and right (For) sides of QD1. In fact, the 
voltages applied to the gates in Fig. la were such that 
rcR ^ Tl- In our model, this is equivalent to having an 
asymmetric t', with t' CI{ ^S> t' L . An analysis of the results 
in this asymmetric regime indicates that the correlation 
between QD1 and QD2 does indeed increase. However, if 
one performs the calculations with the sites in the CR at 
a filling lower than one electron per site (half- filling) , it is 
observed that Si • S2 is gradually suppressed as the elec- 
tron filling falls to a more appropriate level to simulate 
the two-dimensional electron gas in the CR. Although 
one can argue that some of the dependence of the con- 
ductance of QD1 on the charge state of QD2 seen in Fig. 
2 is associated to the correlations between the two dots, 
it is apparent that other effects are also present. This is 
dramatically exemplified by the fact that the cancellation 
of G presented in Fig. 2c occurs for any finite value of t" , 
and of course for t" « 0, one finds that Si • S2 ~ 0. The 
fact that the dip seen in the conductance in Fig. 2c is not 
dominantly caused by correlations between the dots can 
be made more clear by checking the results for the con- 
ductance as Ui (Hubbard interaction in QD1) is reduced 
to zero. In Fig. 5c, results for G are shown for 3 differ- 
ent values of U\ , for the same parameters as for the blue 
curve in Fig. 2c. As U\ decreases from 0.4 (blue) to 0.2 
(green), and then to 0.0 (red), the dip in the conductance 
remains, only becoming narrower, indicating that its ori- 
gin is not associated with many-body interactions, but 
more likely with cancellations typical of T-geometriesi^ 
that occur even in the non- interacting limit. 

In summary, the numerical results qualitatively repro- 
duce the main aspects of important recent experiments 3 
involving nonlocal spin control in nanostructures. The 
main result is that the splitting observed in the ZBA is 
caused by a cancellation in the conductance due to a de- 
structive interference. This so-called Fano anti-resonance 
has its origin in one of the dots being side-connected to 
the current's path. A simple 'circuit model' qualitatively 
reproduces the experiments and offers an alternative to a 
purely RKKY interpretation of the results, underscoring 
that a laboratory realization of the two-impurity Kondo 
system should avoid any geometry susceptible to a Fano 
anti-resonance. 
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